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00 ■ Abstract 

Q : 

Darboux-deformations of short range one-dimensional potentials are constructed by 
means of Gamow-Siegert functions (resonance states). Results include both Her- 
O \ mitian and non-Hermitian short range potentials which are exactly solvable. As 

' illustration, the method is applied to square wells and barriers for which the trans- 

^-H ■ mission coefficient is a superposition of Fock-Breit-Wigner distributions. Resonance 

\ levels are calculated in the long lifetime limit by means of analytical and numerical 

approaches. The new complex potentials behave as an optical device which both 
. refracts and absorbs light waves. 

i : 

g ; 1 Introduction 

I \ Gamow-Siegert functions are solutions of the Schrodinger equation associated to complex 

>• ' eigenvalues and fulfilling purely outgoing conditions [1,2]. Although these functions are 
not square-integrable, they are used to model physical phenomena as resonances which, 

Xf^ . incidentally, are defined experimentally better than theoretically. Some approaches extend 

\ the formalism of quantum theory so that resonant states can be defined in a precise form. 

O I For instance, in the formulation of the rigged Hilbert space, Gamow vectors are studied 

' in the framework of the nuclear spectral theorem [3,4]. Resonances are also connected 

O ■ with the point spectrum of complex-scaled Hamiltonians [5]. 

I In a different context, complex eigenvalues of Hermitian Hamiltonians have been used 

^ I to implement Darboux (supersymmetric) transformations in quantum mechanics [6-9] 

■ (see also the discussion on 'atypical models' in [10]). The transformed Hamiltonians 

include non-Hermitian ones, for which the point spectrum sometimes has a single complex 
eigenvalue [7,8]. This last result, combined with appropriate squeezing operators [11], 
could be in connection with the complex-scaling technique [8]. In general, supersymmetric 
transformations constitute a powerful tool in quantum mechanics [10]. As far as we know, 
however, the study of Darboux-deformations of potentials with resonances is still missing 
in the literature. 

In this paper we deal with one-dimensional short range potentials and their resonances. 
As usual, the complex energy eigenvalue e = E — iV /2 will be taken as a compound of the 
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resonance E and the inverse of the hfetime T. It will be shown that resonances can be 
labelled by a positive integer, according with their closeness to the scattering threshold. 
This integer thus indicates the 'grade of excitation' of the related state: the 'ground 
resonance' will correspond to the least excited Gamow-Siegert state and so on. We shall 
use these solutions to implement Darboux transformations; in this way several exactly 
solvable complex and real potentials will be constructed. As we are going to see, the 
Darboux-Gamow deformation which is produced on the initial potential strongly depends 
on the 'excitation' of the transformation function. Moreover, the complex potentials so 
derived behave as an optical device which both refracts and absorbs light waves. Such 
a behaviour will be illustrated with the simplest short range one-dimensional models: 
square wells and barriers. 

In sections 2 and 3 we present general analytical conditions to construct Gamow- 
Siegert functions for one-dimensional short-range potentials and the corresponding Darboux- 
deformations. In Section 4, we will express the transmission coefficient T as a superpo- 
sition of Fock-Breit-Wigner (Lorentzian) distributions for square well potentials. As we 
shall see, each one of the bell shaped peaks arises since T has a pole at e G C, which 
is a complex eigenvalue of the square well. The Darboux-Gamow deformations of these 
simple wells give rise to complex potentials V = Vii + iVj for which the function Vj shows 
multiple changes of sign, thus each V emits and absorbs flux at the same time. All the new 
potentials of this section have square-integrable functions which represent bound states 
associated with eigenvalues in the point spectrum. In Section 5 the case of square barri- 
ers is analyzed. We shall show that Gamow-Siegert functions also lead to a transmission 
coefficient which is a sum of Fock-Breit-Wigner distributions. In this case, 'haired' square 
barriers appear as a consequence of double Darboux-Gamow deformations for which the 
number of hairs is in close connection with the level of excitation of the transformation 
function. Finally, in Section 6 we finish the paper with some concluding remarks. A short 
appendix is included in order to present further information which, although important, 
can be delayed. 



2 Gamow-Siegert states and one-dimensional short 
range potentials 

Let us consider the Schrodinger equation for a stationary, short range potential V{x), 
defined on the straight line M {h/2m — 1): 

Hu{x, e) = —u"{x, e) -|- V{x)u{x, e) = eu{x, e). (1) 

In the sequel, we shall assume that V{x) is characterized by a strength Vq ^-nd a cutoff 
parameter ( > such that V{x) = if \x\ > (. A bar over z G C stands for complex 
conjugation (z) while zr, zj, denote real and imaginary parts respectively. In general, 
I+{z) denotes the upper half of the 2;-plane {zj > 0) and I-{z) the lower half plane. 
Whenever there is no confusion we use the shortcut notation / = /(x, e, ...), keeping 
implicit the dependence of / on x, e, and other possible variables. 
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It is useful to write the derivative of the logarithm of u as the complex function f3{x, e): 

P:=-^\nu. (2) 
ax 

Let us also write u = $e*", with and S(a;) real functions defined as 

$ = $0 exp ( - / pRdx ) , S = - / Pidx + So (3) 



where $o and So are integration constants included to fix the amplitude $ and phase S 
respectively. Thereby, the current density 

j = i (u'u — uu') (4) 

can be written as 

j = 2S>p = ~2(3i\u\^ = vp. (5) 

Hence, —2/3/ represents the flux velocity f , which plays an important role in the description 
of Gamow-Siegert functions as well as in the Darboux transformations we are going to 
implement. 



2.1 Asymptotic behaviour of the general solution 

The general solution of ([T]) for |x| > C can be written in terms of ingoing and outgoing 
waves: 

u{x < -C) = /e''^" + Le"'^", u{x > C) = iVe"*^" + Se'^"" (6) 

where coefficients /, L, A^, S, depend on the potential parameters C, Vq, and the incoming 
energy fc^ = e. Among these solutions, we are interested in those which are purely outgoing 
waves. The appropriate boundary condition may be written 

lim {u'^iku)= lim {{—(3^ik)u} = {). (7) 

X— >±oo a;— >itoo 

Thus, the second term in each of the functions ([6]) must dominate over the first one. For 
such states, equation (jlj) takes the form: 

i = ±{k + k)\u\'^ , at x^±C. (8) 

If e = -E", then k is either pure imaginary or real according to E is negative or positive. 
If we assume that the potential admits negative energies, we get k± = ±i^y\E\ and (IHl) 
vanishes (the fiux velocity v is zero outside the interaction zone). The solutions 
connected to k+, are not of our immediate interest because they are bounded (indeed, 
the adequate values of /c+ produce square- integrable functions in M). On the other hand, 
antibound states ip^^ ^ increase exponentially as |x| +oo and represent resonances under 
very special conditions (see e.g. [12]). However, these states will not be considered here 
since they do not exist for arbitrary short range potentials. To exhaust the cases of a real 
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eigenvalue e, let us take now E > 0. The outgoing condition ([7]) drops the interference 
term in the density 

p{x;t) = \N\'^ + \S\'^ + 2\NS\cos{2Kx + ATgS/N), x > (, (9) 

so that p = \S\'^ is not an integrable function in neither space nor time (similar expressions 
hold for X < — C)- Remark that flux velocity is not zero outside the interaction zone 
{\v\ = 2\^). Thereby, E > provides outgoing waves at the cost of a net outflow 
j 7^ 0. To get solutions which are more appropriate for this nontrivial j, we shall consider 
complex eigenvalues e. 



2.2 Resonances 

From the discussion of the previous section it follows that, in general, e must be complex. 
We shall write 

e = E-'-T, eR = kl-kj, ej = 2kRkj (10) 
where e = {kn + z/c/)^. According to ([7]), the boundary condition for f3 reads now 

lim {-p ± {ki - ik^)} = (11) 

X— >iboo 

SO that the flux velocity is f+ = +2kFt foTx>C and V- = —2kR for x < Hence, the 
"correct" direction in which the outgoing waves move is given by k^ > 0. In this case, 
the density 

p{x;t) = \u{x,t)\^ = e-^'\u{x)\\ lim p{x;t) = g-^^*""/''*) (12) 

X— >itoo 

can be damped by taking F > 0. Thereby, kj and kji have opposite signs. Since 
kn > has been previously fixed, we have kj < 0. Then, purely outgoing, exponentially 
increasing functions (resonant states) are defined by points in the fourth quadrant of the 
complex fc-plane. As usual, = E will be called a resonance and F = 2|e/| the inverse 
of lifetime. We shall write u{x, e = E — |F) = ipe{x). 

Observe that density f|T2|) increases exponentially for either large |a;| or large negative 
values of t. The usual interpretation is that the compound {ip^, V) represents a decaying 
system which emitted waves in the remote past t — x/v. As it is well known, the long 
lifetime limit (F 0) is useful to avoid some of the complications connected with the 
limit t — oo (see discussions on time asymmetry in [13]). In this context, we shall 
impose the condition: 

r/2 , , 

^ « 1- (13) 

Thus, the level width F is much smaller than the level spacing AE in such a way that 
closer resonances imply narrower widths (longer lifetimes) . An additional condition which 
will be useful in the next calculations reads 

E > AE. (14) 
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In general, the main difficulty is precisely to find the adequate E and T. However, for 
one-dimensional stationary short range potentials and appropriate values of Vq and the 
transmission coefficient T presents a series of peaks. Each sharped peak is a bell-shaped 
curve of width F which can be described by the Fock-Breit-Wigner (FBW) distribution 

-fa-^)- (.„-V4(r/2F 

The center E of the peak defines a wave for which the scattering is ineffective {E is a 
transparency!). Energies which are very close to E, namely \eR — E\ < F/2, induce 
a resonance: the wave interacts with the potential in such a way that the time to cross 
the interaction zone is maximum. The superposition of a denumerable set of FBW distri- 
butions (each one centered at each resonance En,n = 1,2, . . .) entails an approximation 
of the coefficient T such that the larger the number N of close resonances involved, the 
higher the precision of the approximation: 

N 

T ^ UNieR) = ^uj{eR,En). (16) 

n=l 

The above description suggests a graphical method to calculate E and F (see Appendix A). 



3 Darboux-Gamow transformations 

In order to throw further light on the complex function f3 we may note that (j2]) transforms 
the Schrodinger equation ([1]) into a Riccati one 

-(3' + f3^ + e = V. (17) 

Remark that (fT7I) is not invariant under a change in the sign of the superpotential f3 

p' + (3^ + e = V + 2(3'. (18) 

These last equations define a Darboux transformation V = V{x, e) of the initial potential 
V. The Darboux-Gamow deformation is necessarily a complex function 

V = V + 2p' = V-2—\nip,. (19) 
The solutions y = y{x, e, S) of the non-Hermitian Schrodinger equation 

- y" + Vy = Sy (20) 

are easily obtained 

W(^.,^) 

2/oc , 21 

where W(*, *) stands for the Wronskian of the involved functions and ip is eigensolution 
of (lH) with eigenvalue S. 
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Let us analyze the properties of these new solutions. First, if ips is the scattering state 
of H connected with the positive kinetic parameter k, = \fEs-, equations ([6]) become 

V'.< = e'^'^ + L,e-"^\ = Ssc"''' (22) 

where we have assumed a single source at the left of V . Hence, T = \Ss\^ and R = \Ls\'^ 
are such that R + T = 1. The transformed scattering waves are obtained from ( l2Ti) to 
read 

ys< = iP< + if^)e'^^ + {P< - m)L,e-*'^^ = (/?> + m)5,e*'^^ (23) 

Hence, we have 

R=\t\'R, f=\t\X t=^2l±^_ (24) 



After using (fTTj) and considering > we arrive at 



R + T 



l+(lf)' 


1-2^ 




1 + 2^ 



(25) 



where f s 7^ and w 7^ are the flux velocities of the scattering wave ifjg and the Gamow- 
Siegert function ip^ respectively. Thereby, i? + T 1, no matter the value of Vg/v 7^ 1. 
That is, scattering waves ips and their Darboux-Gamow deformations Us share similar 
transmission probabilities. 

Now, let us suppose that Hamiltonian H includes a point spectrum ad{H) C Sp(if). 
If ipn is a (square-integrable) eigenfunction with negative eigenvalue £n, then its Darboux- 
Gamow deformation fl2Tl) is bounded: 



lim = ^{^J En + ik){ lim V'n)- (26) 

Thereby, is a normalizable eigenfunction of if with eigenvalue £n- However, as e is 
complex, although the new functions may be normalizable, they will not form an 

orthogonal set [8] (see also [14] and the 'puzzles' with self orthogonal states [15]). There is 
still another bounded solution to be considered. Function oc (p~^ fulfills equation (!20!) 
for the complex eigenvalue e. Since lim^-^-i-oo |?/eP = e^"^^'^ and fc/ < 0, we have another 
normalizable function to be added to the set {?/„}. 

In summary, one is able to construct non-Hermitian Hamiltonians H for which the 
point spectrum is also ad{H), extended by a single complex eigenvalue Sp(if ) = Sp(if) U 

{4- 



3.1 New short range real potentials 

The iteration of the previous procedure is easy once a set of Gamow-Siegert functions 
{(Pe} belonging to the complex eigenvalues {e} has been given: 

V2 = V + 2(3'^{a) (27) 
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where /?2(«) = P{x,e,a) is a Backlund-like function [16] defined by the finite-difference 
algorithm 

'*<"'^-^<^' + /j(5y^ (^«' 

and fulfilhng 

- (3'2{a)+/3i{a)+a = V, a G C. (29) 
In general V2 = V2{x, e, a) is complex; however, this becomes real if a = e [7,8]: 

''= = ^ + Ki)' = ''"'(v)' P"' 

and the new functions \E' = \E'(x, e, £^), such that if2^ = read 

vl/ = (e-^)^-2(^)y. (31) 

In this case, V2 is a Hermitian short range potential and there is no additional eigen- 
values. Hence, H2 and H are strictly isospectral. A straightforward calculation shows 
that Darboux-Gamow deformations of ipn G L^(]R) correspond to normalizable \E'n while 
scattering states ipg are transformed into \E's, with similar transmission probabilities. 

In the next sections we shall deform the simplest short range one-dimensional poten- 
tials: square wells and barriers. 



4 Square wells 



Let us consider the square well potential 



vix) = -Voe[--\x 



Q{x) 



1, X > 
0, X < 



(32) 



with Vo > and 6 > 0. We shall assume that the only source of particles is at the 
left of the interaction zone (this arbitrariness is allowed by the parity symmetry of the 
potential). In this way, we shall analyze the scattering by an attractive potential of range 
b. 

The scattering solution of the Schrodinger equation ([1]) for potential (|32|) reads: 



u 



^ g-ifc6/2 i^Qg ^ — iq sin ^) sin qx 
+ {q cos ^ — ik sin ^) cos gx] , 



X < --, 



-|<x<| 



(33) 



< X 



where 



kq 



S{k) =[-^] e-^'^" 



(34) 



The kinetic and interaction parameters are defined respectively hj k"^ = e and = k'^ + Vo, 
while A is the function 

A{k) = {k cos f - sin f ) (g cos f - ik sin f ) . (35) 

Let us concentrate on the transmission amplitude S, which is regular in /+ except for 
the points /c+ located on the positive imaginary axis (that is, each fc+ corresponds to a 
bound state, see Appendix A). It is simple to verify that k and —k will produce the same 
peaks in the transmission coefficient, that is, S{k) = S{—k). However, they are in /_ (as 
discussed in Section 2, k is in the fourth quadrant of the fc-plane). Hence, 5* must be 
extended to be analytic in /_. First, let us observe that the relationship 

S{-k)S{k) = 1 (36) 

is fulfilled for the points k which are such that qb = nir. Since —k G /+, this last equation 
defines S{k) = l/S{—k) as an analytical function in /_ but at the points —k which are 
zeros of S in /+. In a similar manner, we realize that S{—k) is analytical in /_ except for 
the points k which are zeros of 5* in /+. In summary, let /c„ be a pole of S in the fourth 
quadrant of the fc-plane, then —kn is also a pole while fc„ and — fc„ are zeros of S. Thereby, 
5* is a meromorphic function of k, with poles restricted to the positive imaginary axis 
(bound states) and the lower half-plane (resonances). On the other hand, if 5* is studied 
as a function of e, it is necessary to consider a two-sheet Riemann surface (fc^ = e) with 
a cut along the positive real axis. Complex poles of S{e) always arise in conjugate pairs 
(corresponding to k and —k) while poles on the negative real axis correspond to either 
bound or antibound states. 



4.1 Analytical approach for calculating resonances 

In this section we shall derive expressions for the resonance E and the width F in terms of 
the potential parameters. We assume that the incoming kinetic parameter is a complex 
number with nontrivial imaginary part. As the presence of the short range potential is 
necessarily encoded in the (complex) interaction parameter, most of the approximations 
will be made on q. 

For complex q, additional to (fTOj) . the real and imaginary parts of e read 

E = ql-Vo- qj, ^ = -2qnqi. (37) 
The combination of equations (137|) leads to: 

\(^iyi'-iVo + Eh-l = (38) 

-2 



where 'y = qj ■ Dropping the negative root we get 



,J . ^'^ f ) (39, 

(V„ + -B) 1 + Jl + i(J^ 



Vo+E 
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On the other hand, since Vq > 0, condition (IT^ produces 



Vo + E E AE 
which, together with condition (fT3|l . imphes 

r/2 

Vo + E 



r/2 r/2 r/2 , , 

< ^ < -r^ (40) 



(41) 



Hence, equation (139|) gives 



This value of qi, when substituted into the expression for r in fl37j) . leads to 

E^qR^- Vo. (43) 

We first notice that the positiveness of E is ensured whenever qn fulfills > Vq. Now, 
a simple comparison between fH3l) and fl371) shows that gf has to be smaller than Vq. We 
then have a relation of order qj < Vq < q\. As we shall see, it is most convenient to 
constrain the values of q^ and qi as weighted by the cutoff 6/2: 

^4)\.<(g)\ HV.,y.^\^. ,44, 

where we have introduced the parameter 9 which encodes not only the range of application 
of our approach but the physical identity of the potential ( |32l) as well. 

Now, making the r.h.s. factor of (!35l) equal to zero, we obtain q = iyT^sin y. After a 
straightforward calculation the corresponding equations for the real and imaginary parts 
are obtained. For 6 » 1 and small values of the last one leads to cos ^ ~ 
Therefore, we arrive at 

^''--^±Uo('] (45) 



2 2 6 



where n is to be determined. After introducing ( l45l) into (1431) . we finally get a discrete 
set of resonances 



E 



n7r\- 
2d) 



71 jr 

K), ^>^^, nen. (46) 



It is clear that n must exceed a minimum value to get a positive energy. Let us take 
n := JT-inf + m, m = 0, 1, . . ., where n^^f is the ceiling function of 26'/7r, i.e., n^^f = j"^] . In 
this way, E^^^ = Eq is the smallest positive energy in fH6|) . Hence 

E^=\''-^^^^i^T-l\V,. »eZ^ (47) 
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The derivation of T is now easy. Making the Lh.s. of A equal to zero we arrive at 
kj = ±a/Vo cos ^ cosh After using (145|) . for small qib/2 one obtains 

kjb 1 / 0/6 \ ^ 

With this last result and fl43l) into (ITU]) , the energy can be rewritten as 

kRbY f Qr^Y , Qi^Y , 1 



2«y V29yn^j^9-' 

the second additive term is smaller than the first one (see equation (jSj)) and 6^^ << 1, 
so equation fH9l) reduces to kji ^ ±\/E. The combined substitution of these last results 
into the expression for F in ffTOj) gives 

2 "5^= IT f^"' 

where the terms which are proportional or greater than the square of qib/2 have been 
dropped (compare our results with those reported in e.g. [17]). 

It is remarkable that F is proportional to f+. For a given value of b {6 fixed), a hasty 
incident particle spreads its peak in T by spending less time in the interaction zone. As a 
result, the peaks centered at hot resonant energies tend to lose the FBW shape by 
excessively overlapping the next ones. The long lifetime limit (F — 0), applied for fixed 
values of 6, is then a good criterion to select velocities as well as potential parameters 
in order to obtain appropriate Gamow-Siegert functions for the square wells. 

4.1.1 Fock-Breit-Wigner distributions 

The above discussion considered a complex kinetic parameter with nontrivial imaginary 
part. As we have seen, approximations leading to E are also useful to get F in the long 
lifetime limit. Next, we are going to verify that E and F so derived correspond to the 
center and width of a FBW peak in the transmission coefficient T. It is convenient to 
rewrite the transmission amplitude (IMl) as follows 

-ikb / 1.2 , 2\ 

Let us consider a small displacement 6E from the resonance energy E^ along the real 
axis of the e-plane. Kinetic and interaction parameters are mapped as — > k + 5k 
and q ^ q + ^SE respectively. Near the resonance we get qb ^ 2nn + ^6E. Thus 
g ^ {2/Tm)SE, where we have used ( 150|) . Hence 



S{k + 5k) 

and finally we get 



e 



-i{k+Sk—n-K)b 



1 - i25E/Tr, 



10 



This last expression is consistent with the previously defined FBW distribution (fT5|) for 
small displacements along the real axis of the e-plane, that is R 9 SE = {er — Em)en^E^- 
A comparison of (H7I) and (l50l) with the results obtained from the graphical method is 
reported in Table [T] of Appendix A. 

4.1.2 Diverse kinds of transformation function 

Let us pay attention to solution fl55]) for a complex kinetic parameter k. Outside the 
interaction zone, it reads 

= e'^^'^e-^'"' + L{k) e-^^«^e'=^^ m> = S{k) e*'=«^e-'=^^ (53) 

According with the discussion of the above sections, 'decaying states' v^e(x) correspond 
to kinetic parameters k living in the fourth quadrant of the complex plane (see Figured]). 
The mirror images in the third quadrant —A;, give rise to 'capture states' ifj behaving 
in the same global way as ip^ (see equations (1531) ). On the other hand, equations (153!) 
evaluated at k correspond to an exponentially decreasing function (see Figure [2]). The 

2 

same is true for —k. We shall write Lp^ = u{x,k ). Although this solution reverts the 
decaying (fT2!) to an exponentially growing dependence on time, it is useful for obtaining 
special cases of Darboux-Gamow deformations, as we are going to see. 




-7.5 7.5 -7.5 7.5 



Figure 1: Real parts of the second (left) and third (right) decaying states of a square potential 
defined by Vq = 50, 6 = 14.2. The Gamow-Siegert function cp^ respectively belongs to the 
complex eigenvalues ei = 3.3029 - i0.5119 and 62 = 6.5823 - iO.7227. 



-10 



10 






-10 





-7.5 7.5 -7.5 7.5 

Figure 2: Imaginary parts of the first (left) and fourth (right) decreasing functions (p-j: of 
the potential described in Figure [TJ The complex kinetic parameters are respectively /cq = 
0.3072 + iO.2484 and A?3 = 3.1496 + i0.2811. 
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4.2 New complex and real wells 

It is common to use complex potentials V = Vr + iVi, Vj = const, for describing either 
absorption or emission of flux according with the sign of Vj (see e.g. [18]). In contrast, the 
case when Vj is an arbitrary function of the position is rarely analyzed because the related 
eigenvalue equation is involved. However, in Section 3 we realized that Darboux-Gamow 
deformations V are non-trivial exactly solvable complex potentials. 




Figure 3: The real (left) and imaginary (right) parts of the Darboux-Gamow deformation V of 
a square well (dashed) for which Vq = 16 and 6 = 5. The lowest decreasing function ip^ of the 
initial potential is the transformation function with ko = 1.7504 -|- iO.7657. 

A typical Darboux-Gamow deformation of the square well (132|) is shown in Figure [31 
First observe that the initial parity symmetry is now broken. It is also remarkable the 
presence of several local minima and maxima as well as multiple changes of sign in both 
functions Vr and Vj. In a simple optical model the incident beam crosses a sequence of 
'obstacles' (potential's maxima and minima), each one absorbing a fraction of the beam 
and performing some operations on the rest in such a way that the states are never mixed 
(see discussions on the optical bench in [10]). As a consequence, the 'optical' potential 
V is no longer self-adjoint and corresponds to an optical device which both refracts and 
absorbs light waves. Moreover, this non-Hermitian potential has a point spectrum and 
normalized eigenfunctions which inherit its broken parity symmetry (see Figure H]). 



0.5 




-2.5 2.5 -2.5 2.5 



Figure 4: Left: The probability density belonging to the normalized ground state (dashed) of 
a square well potential (Vq = 16, 6 = 5, see Figure [8]) and its Darboux-Gamow deformation 
(continuous curve). The transformation parameters are the same as in Figure El Right: The 
related second excited states. 

On the other hand, new real short range potentials are obtained by simultaneously 
using (f^ and (p-^ in (p71) - (!28|) . The result is illustrated in Figure [51 with a real potential V2 
for which Vq = 16 and 6 = 5. It is notable that local deformations (that is, the number of 
maxima and minima) increases with the level of excitation of the transformation function. 
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-2.5 2.5 



— 













-2.5 2.5 




Figure 5: Twice Darboux-Gamow deformed square well potential (Vq = 16, 6 = 5) for the same 
transformation parameters as in Figure [3l A detail of the bottom of this potential is shown at 
the right. 



Finally, let us remark that, although our one- dimensional model is very simple, exten- 
sions to s- waves in three dimensions are immediate by properly selecting the odd solutions 
and adding to the initial potential an impenetrable wall along the negative axis. 



5 Square barriers 



The barrier potential (Vq > 0) 



V{x) 



\x 



(54) 



is usually studied clS db special case of ([32 



a negative number. Thus, kinetic and interaction parameters now read k"^ 
— k'^ — Vo- The involved solutions are then obtained by taking Vq —>■ —Vq and e 



for which the strength Vq is allowed to be 

e and 

= fc' — Vq- i Jie mvolved solutions are tlien obtamed by taking Vq —>■ — Vq and e = E > 
in fl33l) - fl35l) . However, a simple inspection of (H7j) and fl50l) shows that this is not the case 
in analyzing resonances {E^n and F^ are respectively negative and pure imaginary numbers 
if Vq is simply negative!). In this section we shall obtain the appropriate expressions for 
the complex energies flTU]) . As in the previous calculations, we look for the complex roots 
of A = fulfilling conditions (fT3l) and (fT^ . A straightforward calculation on the l.h.s. of 
( |35|) for the new k and q leads to 



, 777 . QRb . . Qib 
g/j = ±VVoSin — smh— , qi 

The smallness of the weighted qj is now ensured by taking 

QRb {2s + l)n 



: /TT iBb qib 

± V I/q cos cosh 



5. 



0,1,2,... 



(55) 



(56) 



where 5 is a slight perturbation of the weighted q^. The introduction of (!55|) and (156!) 
into the real and imaginary parts of e (compare with equation (!37j) ) gives 



E„ 



Vo 



(2s + l)7r 



+ 1 



2 



9 



0,1,2, 



(57) 
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In a similar form, the r.h.s. of equation (1351) leads to expressions for E and F which are 
the same as (!57|l but they hold for semi-integer values of s. Thus we can write En and 
r„ with n = 1,2,3,... such that s = {n — l)/2. Therefore, En > Vq for any integer 
n > 1. Moreover, is always a positive number which grows as n — oo. Thus, resonant 
behaviour will be better distinguished, and it will have longer lifetime, for energies which 
are close to the top of the barrier. In other words, the resonance threshold is established 
by the strength of the potential (compare our results with those reported in e.g. [19]). In 
Table [2] of Appendix A we report a comparison of the results (|57|) with those obtained by 
means of the graphical method. 



10000 
5000 


-5000 
-10000 




-15000 -5000 2500 



Figure 6: The Argand-Wessel diagram of the Darboux-Gamow deformed square barrier for 
Vq = 1000, 6 = 5 and its first resonance (see Table [2]). The external rings of this cardiod-like 
curve correspond to large negative values of the position x. 

Figure [6] shows the Argand-Wessel diagram of the Darboux-Gamow deformation of a 
square barrier. The real and imaginary parts of this new potential show a series of maxima 
and minima which, in a similar manner as for the complex wells, can be modelled by the 
optical bench described in [10]. Once the method is iterated by means of ip-^, the twice 
Darboux-Gamow deformed barriers are real. These new barriers present 'hair' over the 
top which induces stronger resonant phenomena (see Figure [7]). Let us remark that the 
number of hairs increases with the excitation of the transformation function, just as it 
occurs for the nodes in bounded wavef unctions. Thus, the lowest resonance induces only 
one very localized distortion (a couple of hairs combed in opposite direction), the second 
one induces four hairs and so on. Observe also the global asymmetry of the new barriers, 
that is, V2{—x) 7^ V2{x). Finally, it is reasonable to assume that these haired barriers 
will induce delays on the scattering states which are longer than the delay produced by 
conventional square barriers (see simple models in [20]). The analysis of this phenomenon 
win be published elsewhere. 

6 Concluding remarks 

We have derived analytical expressions for the resonant energies of short range one- 
dimensional potentials in the long lifetime limit of the resonance levels. The same calcu- 
lations leading to the 'quantization' Em of the resonances are appropriate to obtain the 
corresponding lifetime inverses F^. These results allowed the description of the transmis- 
sion coefficient as a superposition of Fock-Breit-Wigner distributions, which is as good as 
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-4 4 

Figure 7: The top of twice Darboux-Gamow deformations of the square barrier Vq = 1000, 
b = 10. The lowest decreasing function produces a very locahzed distortion (couple of 'hairs') 
on the new top (gray curve). The second (dotted curve) and third (sohd curve) ones produce 
respectively two and three distortions (see Table [2] for specific values of the complex eigenvalues). 

large is the number of resonances involved. Moreover, they also were used to transform 
Hermitian short range potentials into non-Hermitian ones which either preserve the ini- 
tial energy spectra or include a complex eigenvalue with square-integrable eigenfunction. 
The new complex potentials are 'opaque' in the sense that they simultaneously emit and 
absorbs flux, just as an optical device which both refracts and absorbs light waves. The 
iteration of the method (e and e combined) produces new real potentials for which the 
energy spectrum is exactly the same as the initial one. 

The Darboux-Gamow deformations of square wells and square barriers have been 
presented as straightforward applications. In the firs case, it has been shown that square- 
integrable eigenfunctions are mapped into square-integrable deformed eigenf unctions. In 
the second case we obtained 'haired' square barriers as a consequence of double Darboux- 
Gamow transformations. The hairs are very localized deformations on the top of the new 
barriers and are such that their number depends directly on the excitation of the involved 
transformation function. In all cases, the scattering states preserve their global properties 
after the deformations. 

Some final comments would be important. Processes in which the incident wave falls 
upon a single scatterer are fundamental in the study of more involved interactions [21]. 
In general, for a single target the scattering amplitude is a function of two variables (e.g. 
energy and angular momentum). Our model corresponds to the situation in which one 
of the variables is held fixed (namely, the angular momentum). A more realistic three 
dimensional model is easily obtained from our results: even functions are dropped while 
an infinitely extended, impenetrable wall is added at the negative part of the straight line. 
Such a situation corresponds to s-waves interacting with a single, spherically symmetric, 
square scatterer (see e.g. [18]). On the other hand, the production of confined states 
within the continuum associated with lattice impurity (contact effect) has been recently 
reported in the context of Darboux-deformations [22]. Such a result indicates that our 
Darboux-Gamow approach can be extended to the case of spatially oscillating square well 
potentials. Then, the presence of 'locahzed' lattice deformations produced by resonant 
states could be in connection with the isolated transitions, observed in semiconductors, 
from a bound state within a quantum well to a bound state at an energy greater than the 
barrier height [23]. Work in this direction is in progress. 
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Appendix A. 



T4) = 1000, 6 = 20, ninf = 202, e = £;,„- ir„/2 


Graphic 


Analytical 


06.798680-^0.527888 


06.798344 -i0.521472 


TO = 


16.790717 -i0.843441 


16.791319-^0.819544 


TO = 1 


26.832718-^1.084738 


26.833641 -^1.036023 


TO = 2 


36.926378-^1.295156 


36.925312-^1.215324 


TO = 3 


47.065373-^1.489006 


47.066330-^1.372098 


TO = 4 


57.258014-^1.673245 


57.256697-^1.513363 


TO = 5 


67.497189-^1.852008 


67.496412-^1.643124 


TO = 6 


77.786734-^2.028109 


77.785474-^1.763421 


TO = 7 




^0 = 743, 6 = 22, ninf = 191, e = S„ - ir„,/2 


Graphic 


Analytical 


00.911072-^0.174723 


00.911235-^0.173561 


TO = 


08.720839-^0.547966 


08.721274 -i0.536941 


TO = 1 


16.571405-^0.768964 


16.572095-^0.740160 


TO = 2 


24.463107-^0.951658 


24.463700-^0.899287 


TO = 3 


32.397141-^1.116079 


32.396089-^1.034864 


TO = 4 


40.370347 -il. 270414 


40.369261 -^1.155214 


TO = 5 


48.381943 -il.419086 


48.383217-^1.264691 


TO = 6 


56.436735-^1.564930 


56.437956-^1.365912 


TO = 7 



Table 1: Representative numerical results for the resonance energies and widths of a square well. The 
potential parameters have been chosen such that the function ujNi^R.) matches the transmission coefficient 
T, at least for the first 8 peaks. 

The physical properties of the solution fl33l) are encoded in A. For example, the 
product uA makes evident that, after suppressing the incoming wave (i.e., A = 0), still 
there are 'reflected' and 'transmitted' waves as well as a finite wave inside the interaction 
zone. Thus, the outgoing condition ([71) is automatically satisfied. Moreover, since specific 
values of A;+ = i>y\E\ produce both the quantization of E' < and the vanishing of 
the related solutions at ±oo, the usual interpretation of bound states as resonances with 
infinite lifetime (F = 0) is then recovered. Conventional expressions for these solutions 
are usually chosen to be simultaneously eigenstates of the parity operator, i.e. odd and 
even functions: 

x<-| 

l<x<l (58) 



V'odd(a;) = < 



exp{^y\E\ x), 

sin| ) ^1^^^' 

— exp(— a/|-E| x). 
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Vo = 1000, 6 = 5, e^En- ir„/2 


1 T" o It 1 
Vjrl cLjJlllC 


Analytical 




1000.394784- i0.009987 


n — 1 


1001.579136-^0.040036 


1001.579136 -i0.039949 


n = 2 


1003.553057-^0.090291 


1003.553057-^0.089886 


71 = 3 


1006.316546-^0.161043 


1006.316546-^0.159797 


rt = 4 


1009.869604-^0.252698 


1009.869604-^0.249683 


n = 5 


1014.212230-^0.365794 


1014.212230-^0.359544 


n = Q 


1019.344424 -iO. 501014 


1019.344424-^0.489379 


n = 7 


1025.266187-^0.659204 


1025.266187-^0.639189 


71 = 8 




K) = 1000, 6 = 10, e = S„ - ir„/2 


Graphic 


Analytical 


1000.098696-^0.001248 


1000.098696-^0.001248 


71 = 1 


1000.394784 -i0.004996 


1000.394784-^0.004993 


n = 2 


1000.888264 -i0.011248 


1000.888264-^0.011235 


n = 3 


1001.579136-^0.020013 


1001.579136-^0.019974 


= 4 


1002.467401 - i0.031303 


1002.467401-^0.031210 


ft = 5 


1003.553057 -i0.045134 


1003.553057-^0.044943 


?i = 6 


1004.836106 -i0.061525 


1004.836106-^0.061172 


n — 7 


1006.316546-^0.080501 


1006.316546-^0.079898 


71 = 8 



Table 2: Representative numerical results for the resonance energies and widths of a square barrier. 
The potential parameters have been chosen such that the function (^^{eii) matches well the transmission 
coefficient T, at least for the first 8 peaks. 



exp{^y\E\x), X < — 2 



cxp(-f Vl^l) \ „„„ „^ b ^ ^ ^ b 



COS -2 



jcosqx, -|<a;<| (59) 



V'even 



2 

The energy eigenvalue conditions for these solutions are respectively 



exp{—^y\E\ x), I < X 



,2 



cotg=—— , ta.ng = — (60) 

g g 

with 6 defined in (jUj) and g = qb/2. Figure [8] shows some of the first normalized eigen- 
functions we are dealing with. 

On the other hand, for one-dimensional short range potentials the transmission ampli- 
tude S could be interpreted as a 1 x 1 scattering matrix. So, to measure the cross section 
as a function of the energy of incoming particles corresponds to identify the local maxima 
of the transmission coefficient T. First, notice that a sharped peak of T can be connected 
with a Fock-Breit-Wigner distribution ( |T5i) . just as it is illustrated in Figure [9l The peaks 
surrounded by a couple of local minima which are above 1/2 will be dropped. The top of 
each of the included peaks, projected on the energy axis, defines a resonance. The width 
r of each peak is then defined by the distance between its right and left intersections with 
the horizontal at 1/2. Function ujN{eR) in (fT6l) includes N of these peaks, counted from 
left to right. In Figure [9] it is also shown a case for which the matching between T and 



17 




-0.5 \ _/ y 

-2.5 2.5 

Figure 8: The ground state (continuous curve) and the first two excited bound states of a 
square well Vq = 16, 6 = 5. The energies are respectively = —15.6379, E\ = —14.6983, 
E2 = -13.0812. 

cujv is very good for the first ten resonances. Finally, Table [1] shows the results obtained 
with the graphical method as compared with the analytical procedure of Section 4 for 
the square well and different values of the range b and strength Vq. Table [2] shows the 
corresponding results for the square barrier of Section 5. 




2 4 6 a 10 1000 1002 1004 1006 1008 1010 



Figure 9: Left: First resonance of the square well Vq = 992.25, b = 20. The graphic definition 
of the width F corresponds to the distance between the vertical lines centering the peak. Right: 
Functions T and con (dotted curve) for the square barrier Vq = 10^, b = 10. In this case, the 
FBW sum matches well the transmission coefficient for the first ten resonances (see Table [2]) . 
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